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Abstract 

In this article we analyze the nuclear matrix elements (NME) of the neutrinoless 
double beta decays of the nuclei ^^Ca, ^^Ge, ^^Se, i24sn, ^"^^Te, ^^^Te and 

^^^Xe in the framework of the Interacting Shell Model (ISM). We study the relative 
value of the different contributions to the NME's, such as higher order terms in the 
nuclear current, finite nuclear size effects and short range correlations, as well as 
their evolution with the maximum seniority permitted in the wave functions. We 
discuss also the build-up of the NME's as a function of the distance between the 
decaying neutrons. We calculate the decays to 0^^ final states and find that these 
decays are at least 25 times more suppressed with respect to the ground state to 
ground state transition. 

Key words: Shell Model, Double beta decay matrix elements. 
PAGS: 23.40.Hc, 21.60.Cs, 27.40.+Z, 27.50.+e, 27.60. +j 



1 Introduction 



The discovery of neutrino oscillations in recent experiments at Super-Kamiokande 
[1], SNO [2] and KamLAND [3] has changed the old conception of neutrinos 
by proving that they are massive particles. According to the origin of their 
mass, neutrinos can be either Dirac or Majorana particles, the latter case be- 
ing particularly interesting since it would imply an extension to the standard 
model of electroweak interactions. In this scenario neutrinos are their own an- 
tiparticles and lepton number is not conserved. Besides, it happens that the 
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best way to detect one of these violating processes and consequently to estab- 
lish the Majorana character of the neutrinos would be the observation of the 
neutrinoless double beta decay (Oz//?/3). 

Double beta decay is a very slow weak process. It takes place between two 
even-even isobars when the single beta decay is energetically forbidden or 
hindered by large spin difference. Two neutrinos beta decay is a second order 
weak process — the reason of its low rate — , and has been measured in a few 
nuclei. The Ou^jS decay is analog but needs neutrinos to be Majorana particles. 
With the exception of one unconfirmed claim [4f5] . it has never been observed, 
and currently there is a number of experiments either taking place [6f7f8] or 
expected for the near future — see e.g. ref. [9] — devoted to detect it and to 
set up firmly the nature of neutrinos. 

Furthermore, Oi>f3(3 decay is also sensitive to the absolute scale of neutrino 
mass, and hence the mass hierarchy — at present, only the difference between 
different mass eigenstates is known. Since the half-life of the decay is de- 
termined, together with the neutrino masses, by the nuclear matrix element 
(NME) for this process, the knowledge of these NME's is essential to pre- 
dict the most favorable decays and, once detection is achieved, to settle the 
neutrino mass scale and hierarchy. 

Two different and complementary methods are mainly used to calculate NME's 
for Oz//3/3 decays. One is the family of the quasiparticle random-phase approxi- 
mation (QRPA). This method has been used by different groups and a variety 
of techniques is employed, with results for most of the possible emitters [TO|TI] . 
This work concerns to the alternative, the interacting shell model (ISM) [T2] . 

In previous ISM works [T3lfT4] . the NME's for the Ou^P decay were calculated 
taking into account only the dominant terms of the nucleon current. However, 
in ref. [I5] it was noted that the higher order contributions to the current 
(HOC) are not negligible and it was claimed that they could reduce up to 
20%-30% the final NME's. Subsequently, other QRPA calculations [TMT] have 
also taken into account these terms, although resulting in a somewhat smaller 
correction. These additional nucleon current contributions have been recently 
included for the first time in the ISM framework |18) . and it is the aim of this 
work to look at them in more detail in order to grasp their relevance for the 
NME's. Then we give a comprehensive overview of the present state-of-the-art 
ISM and QRPA results by comparing the relative contribution of the different 
contributions to the NME's. In particular, the short range correlations (SRC) 
are modeled both by the Jastrow prescription and by the UCOM method 
[To] . The radial dependence of the NME — first discussed in ref. [20] — is also 
studied within the ISM. 

In addition to the usual calculation of Oz//5/3 decays to the ground states of the 
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final nuclei we have also computed the corresponding NME's for the decays 
to the first excited 0^ states. This comparison has already been performed 
for most emitters using QRPA methods |21f22f23f24f25| . concluding that the 
transitions to excited states are much more suppressed. However, it is inter- 
esting to test these results from the ISM point of view, and even to improve 
them, since in refs. [2T|22|23I[24] the new hadronic current contributions were 
neglected and in ref. |25| the short range correlations were not treated properly 
(see ref. [26] )• 



2 Theoretical framework 

The starting point for the Oz//3/3 decay is the weak Hamiltonian: 

Hw = -^ {jL.J^^) + h.c, (1) 



where is the leptonic current, and the hadronic — nuclear — counterpart 
is given in the impulse approximation by: 



9v (g^) r - ^9M (g^) 1^ 



-^7A(g')7^75 + ^7p(g')g^5)^, (2) 



with q'^ the momentum transferred from hadrons to leptons, this is, = 

Pneutron Pproton- 

In the non relativistic case, and discarding energy transfers between nucleons, 
we have: 

Ji' (x) = E r- {9'' J' {q') + 9''Jl {q')) S (x - r„) , (3) 

n=l 



where: 

J' {q')=9v {q'), 

Jn[q)= 19 M [q ) + 9A[q) a„, -9p[q) (4) 
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The parametrization of the couplings by the standard dipole form factor — 
to take into account the finite nuclear size (FNS) — and the use of the CVC 
and PCAC hypotheses — for the magnetic and pseudoscalar couplings Qm and 
gp — are those described in ref. [I5]. We take as values of the bare couplings 
gv (0) = 1 and gA (0) = 1.25. 

Due to the high momentum of the virtual neutrino in the nucleus — «ilOO 
MeV — we can replace the intermediate state energy by an average value and 
then use the closure relation to sum over all the intermediate states. This 
approximation is correct to better than 90% [27] . 



We also limit our study to transitions to 0"*" final states, and assume electrons 
to be emitted in s wave. Corrections to these approximations are of the order 
of 1% at most, due to the fact that in the other cases effective nuclear operators 
of higher orders are needed to couple the initial and final states. 

With these considerations, the expression for the half-life of the Oz//?/3 decay 
can be written as I28ll29l: 



(T°;f (0+ ^ 0+))"' = Goi \m''^^^ f (5) 



where {m^) = \ J2k U^k^kl is the effective Majorana neutrino mass, a combina- 
tion of the neutrino masses due to the neutrino mixing matrix U and Gqi 
is a kinematic factor — dependent on the charge, mass and available energy of 
the process. M^^f^^ is the NME object of study in this work. As we see, the 
neutrino mass scale is directly related to the decay rate. 

The kinematic factor Gqi depends on the value of the coupling constant gA- 
Therefore the NME's obtained with different gA values cannot be directly 
compared. If we redefine the NME as: 



These M'^^^^'s are directly comparable no matter which was the value of 
gA employed in their calculation, since they share a common Gqi factor — 
that of gA = 1.25. In this sense, the translation of M'^'^^^'s into half-lives is 
transparent. The QRPA results obtained with different gA values are already 
expressed in this way by the authors of refs. [20f26] while the results of refs. 
[T6|17II30| are not. The latter have been translated into the form of eq. ([6]) 
when comparing to our results. 
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Looking back at the NME, it is obtained from the effective transition operator 
resulting of the product of the nuclear currents: 

O (q) = -h^ (g) + /^GT ^^^^ _ s^^, (7) 



where S^^ = 3 (qcTnqcrm) — (Jn^'m is the tensor operator. The functions h (q) 
can be labeled according to the current terms from which they come from: 



h^"" (?) = h^I (?) + C (?) + (?) + hZ. (?) 
(?) = hip (g) + h^p (g) + h^^ (g) , 



whose explicit form can be found in ref. [15] . 

In ISM works previous to ref. [18], only haa and hyy terms were considered. 
However, rough estimates of the value of all terms taking q ~ 100 MeV give: 
haa ~ hjjy ^ 1, hap ~ 0.20, hpp ~ 0.04 and hmm ~ 0.02. Therefore, according 
to this figures, certainly hap cannot be neglected. Since the Gamow- Teller 
contribution will be the dominant one, and both the hpp and hmm have the 
same sign and opposite to hap, it seems sensible to keep all these terms in the 
calculation. In Section [3] we will see how these terms indeed contribute to the 
total NME. 

Integrating over q we get the corresponding operators in position space, which 
are called the neutrino potentials. Before radial integration they look like: 



2 R 
2 R 





oo 



h^/^^ (q) 
Jo(?r) f , \ qdq, 
{.l + N 

-32 (qr) ^^^^^qdq, 



(9) 



where j„ (x) are the spherical Bessel functions, r is the distance between nu- 
cleons and R, which makes the result dimensionless, is taken as R=rQA^/^, 
with ro = 1.2 fm. 

Finally, the NME reads: 
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' gy (0) ' 
.9 A (0), 

I E r„ {-V^ (r) + (r) a^a^ - (r) Sj; J |o+) (10) 



In the calculation of the NME we need to take into account the short range 
correlations. Using a standard Jastrow type function, the NME's are modified 
as I3T1I32] : 



{0^\V (r) |0+).., = (0+/ (r) \V (r) \f (r) 0+) 

= {0^\firfVir)\Ot), (11) 

with / (r) = 1 — e~"^^ (1 — 6r^), where a = 1.1 fm~^ and b = 0.68 fm~^. 



3 The NME's disassembled 



Following these considerations, we have performed calculations for the 0^(3(3 
decays of the emitters ^^Ca, ^^Ge, ^^Se, ^^^Sn, ^^sje, ^^^Te and ^^^Xe, using 
the ISM coupled code described in ref. [T2j, ideally adapted for the calculation 
of 0"*" states. Full diagonalizations are accomplished within different valence 
spaces and effective interactions. For instance, the decay of ^^Ca is studied 
in the pf major shell, and the KB3 interaction is employed. For the case of 
^^Ge and *^Se, the valence space consisting on lp3/2, 0f5/2, lpi/2 and 0g9/2 is 
diagonalized using the GCN28.50 interaction. Finally the 0g7/2, ld3/2, ld5/2, 
2si/2 and Ohii/2 valence space and the GCN50.82 interaction are used in the 
decays of ^^"^Sn, ^^'^Te, ^^°Te and ^^^Xe. The latter spaces and interactions will 
be discussed in detail elsewhere [33] . 

QRPA valence spaces comprise at least two major oscillator shells. The effect 
of the orbits excluded in our ISM calculations was studied in Ref. [34], in 
the particular cases of A = 82 and A = 136. Proton and neutron 2p-2h 
excitations were considered separately, obtaining an increase in the NME's 
not larger than 20%. However, this number may be an overestimation since 
the occupancies obtained for the extra orbits are larger than that of QRPA 
calculations. Furthermore, when more than one orbit was added to the valence 
space, the NME increased in the same amount as with only one extra orbit, 
showing that different contributions do not sum up. An extended study of 
these effects is currently in preparation. 

The results of the Shell Model calculations were shown in ref. [H], where a 
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Table 1 

Evolution of the ^^Ge — > ''^Se NME's described in Section [2] as we successively add 
the ap, pp and mm HOC contributions. 





vv + aa 


+ap 


+PP 


+mm 


Mhoc 


4.04 


2.82 


3.29 


3.29 


Mhoc+fns 


3.45 


2.49 


2.80 


2.96 


■Mfioc+src 


2.85 


2.12 


2.36 


2.36 


M0UI3I3 


2.70 


2.01 


2.21 


2.30 






-25% 


+7% 


+3% 



comparison of the final values to the ones from recent QRPA calculations v^^as 
made, pointing out the importance of including all the seniority components 
in a given space. Here v^^e will have a closer look at the NME's, in order to gain 
insight into the details of the calculation and recognise possible uncertainties. 



In Tabled] we show the contribution of the additional nucleon current terms in 
the particular case of ^^Ge — ^ ™Se OuPP decay. In addition to the final result we 
also give the partial ones not considering FNS and/or SRC. We see that these 
terms contribute to reduce the bare matrix element in approximately 20%. 
However, as was pointed out in ref. [H], the effect of the HOC contributions 
in the full NME — FNS and SRC taken into account — is a bit smaller, 15%, 
which only means that HOC contributions are slightly more regularized by 
FNS and SRC than the lower order — aa and vv — terms. In addition to that, 
the mm contribution to HOC, which increases the NME and hence reduces 
the net HOC effect, vanishes in the bare case. The behaviour of all other 
decays is very similar. From this table we also conclude that the overall HOC 
contribution is of the expected relative value, and, moreover, the individual 
terms are also very close to their estimate in Section O they are found only a 
bit larger, but with the correct relative figures. 

It has recently been discussed in ref. [18] that the pairing interaction favours 
the OuPP decay and that, consequently, truncations in seniority tend to over- 
estimate the value of the NME's. This is, when we perform truncations in 
seniority in the initial and final wave functions, the value of the full NME 
is not converged, in general, until the full seniority calculation is done. This 
may explain the difference of the results — due to their different treatment 
of pairing type correlations — for the NME's obtained between the ISM and 
QRPA approaches. 

High seniority components are strongly connected to nuclear deformation. As 
an example, we show in Table [2] the decomposition of the wave function of the 
nucleus ^^Ge — that would exhibit a ficticious OuPP decay to its mirror ^^Se — 
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Table 2 

Decomposition of the wave function of the ground state of ^^Ge according to its 
seniority components, in percentage, for different values of the deformation /?. 







s = 


s 


= 4 


s = 6 


s = 8 


s = 10 


0.15 




78 




20 


1 


1 





0.20 




39 




43 


7 


10 


1 


0.25 




20 




43 


14 


20 


3 


0.30 




6 




32 


21 


31 


10 


Table 3 

Comparison of the NME for the ^"^Ge '^^Se Oi//3/3 decay for this work (ISM) and the 
Jyvaskyla (JY07) and Tubingen (TU99) groups. The values of TU99 were originally 
calculated with vq = 1.1 fm, and have been corrected to be directly comparable with 
the others. 




-^bare 


Mfns 


Mhoc 


Mh.+fns M^'^CP 


C/ b/fns 
/('hoc 


%fns %src 


ISM 


4.04 


3.45 


3.29 


2.96 2.30 


19/14 


10 22 


JY07 p] 


8.53 




7.72 


6.36 4.72 


9/- 


18 26 


TU99 [I5] 




7.03 




0.63 


-/20 


^10 ^20 



for different deformations, obtained by adding a variable extra quadrupole- 
quadrupole term to the interaction. We see that, as the nucleus becomes more 
deformed, the high seniority components become more important. 

We can check if the discrepancy of the NME's remains in the HOC, FNS and 
SRC contributions. For that purpose, in Table [3] we compare in detail the 
NME's for A = 76 for ISM and QRPA calculations, enclosing the partial re- 
sults as well. We also append the amount of correction due to each HOC, FNS 
and SRC to the previous step in the calculation. The SRC are always of the 
Jastrow type unless otherwise stated; we will return to this point later. Similar 
figures are found for any other decay. If we focus in the relative importance 
of the different approximations, we can see that all three calculations show 
the same trend, though the details may change a bit from one to another. 
Hence we can conclude that HOC contributions reduce the NME about 10- 
20%, taking into account FNS effects produces an additional 10-20% decrease 
and finally Jastrow type SRC reduce the NME 20-25%. The overall effect from 
the original bare NME adds up to 35-45%. In our calculations, these three con- 
tributions when applied to the bare NME, give reductions of 20%, 15% and 
30%, respectively — Mgrc = 2.85, not included in Table [3l Therefore, these 
effects do not pile up, the total reduction amounting to roughly 70% of the 
sum of the individual contributions. 
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Table 4 

Evolution of the NME with the maximum seniority (s^) permitted in the wave 
functions of ^^Ge and ^^Se, including different contributions to the full operator. 



Sm 




Mfns 


Mhoc 


Mh, -(_ f r) s 

It* I /to 






J ""^ 


9o g'pQ 





12.31 


11.16 


10.49 


9.83 


8.59 


1 5 


''12 


6 


13 


4 


8.84 


7.87 


7.44 


6.89 


5.82 


16/ 


'12 


7 


16 


6 


8.01 


7.11 


6.73 


6.22 


5.23 


16/ 


'U 


8 


16 


8 


5.63 


4.90 


4.66 


4.25 


3.34 


17/ 


'12, 


9 


21 


10 


4.64 


4.00 


3.81 


3.45 


2.74 


18/ 


'U 


9 


21 


12 


4.10 


3.50 


3.34 


3.01 


2.34 


19/ 


'U 


10 


22 


14 


4.04 


3.45 


3.29 


2.96 


2.30 


19/ 


'lA 


10 


22 



Table 5 

Evolution of — as a function of the maximum seniority allowed in the wave 
functions, Sm, for all the studied decays. 





A = 48 


A = 76 


A = 82 


A = 124 


A = 128 


A = 130 


A = 136 





0.33 


0.31 


0.30 


0.27 


0.27 


0.27 


0.26 


4 


0.19 


0.23 


0.21 


0.15 


0.20 


0.19 


0.15 


6 


0.17 


0.22 


0.20 


0.16 


0.20 


0.19 


0.16 


8 


0.16 


0.16 


0.13 


0.15 


0.16 


0.15 


0.15 


10 




0.14 


0.12 


0.15 


0.15 


0.15 




12 




0.12 


0.11 


0.15 


0.15 






14 




0.12 













Therefore, contrary to what happens for the NME's, we see that the relative 
contributions of HOC, FNS or SRC are similar for ISM and QRPA calcu- 
lations, which seems to point out that their contribution to the full NME is 
not affected by seniority truncations. This is confirmed in Tabled! where these 
partial contributions as a function of the seniority are shown. The ^^Ge — > ''^Se 
Oz//5/3 decay is chosen again, but the same conclusion is obtained for all other 
transitions. Even though a small decrease of all the contributions is seen, at 
the = 4 level — leading order in QRPA — the relative value of HOC, FNS 
and SRC is essentially that of the full calculation. 

Another aspect of the NME that could be sensitive to the treatment of the 
pairing correlations is the ratio of Fermi to Gamow- Teller terms by means 

of the coefficient = (f^lsy) M^/M^'^, represented in Table [5] for all the 
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Table 6 

Comparison of the values of of this work (ISM), this work with seniority = 4 
and the QRPA results of the Jyvaskyla (JY07) and Tiibingen (TU07) groups. The 
TU07 result is taken prior to SRC correlations, see ref. [26] , 







ISM 


ISM {Sm = 4) 


JY07 pfTT] 


TU07 [15] 


™Ge 




0.12 


0.23 


0.27 


0.32 




82Kr 


0.11 


0.21 


0.26 




i28Te 


i30xe 


0.15 


0.20 


0.31 




i30Te 




0.15 


0.19 


0.31 


0.36 


i36Xe 




0.15 


0.15 


0.27 





studied nuclei as a function of the seniority. Unlike the precedent case, here 
we see that correlations affect in a different manner to these contributions, in 
such a way that the ratio decreases as we allow higher seniority components 
in the wave functions. This trend is not seen, however, in the A = 48, A = 
124 and A = 136 cases. But these are precisely the nuclei for which the low 
seniority truncation works better [18], being emitters which, in their natural 
valence spaces, only consist on neutrons ("^^Ca and ^^^Sn) or protons (^^^Xe), 
leading to wave functions dominated by low seniority components. We compare 
our full results and those truncated in seniority with the QRPA figures in Table 
El observing that ISM values are smaller than QRPA's, but the truncated 
= 4 results are always closer to them. This is to say, one may attribute 
the discrepancy in x^ between ISM and QRPA to the seniority truncations, 
as was the case for the complete NME. 

So far we have argued that there is agreement in the relative importance of 
the different pieces of the Oz//?/? NME studied, i.e. HOC, FNS and SRC —of 
Jastrow type — between the ISM and QRPA in its different versions, while 
the difference in x^ be attributed, as in the case of the full NME, to 
the fact that high seniority components may not be fully included within the 
QRPA approach. However, there are still two issues regarding which there are 
discrepancies between the different methods and authors. 

On one hand there is the role of tensor part of the NME, which we can quantify 
similarly to the Fermi case by the ratio to the Gamow- Teller contribution, 
= M^/M*^-^. While for some QRPA authors [H] this quantity amounts 
to 5% for A = 76 and 8% for A = 130 —prior to SRC, after that their 
influence is necessarily larger since tensor type contributions, contrary to Fermi 
and Gamow- Teller parts, are not reduced by SRC — , others claim that their 
contribution is negligible [16]. Our ISM results are shown in Table [7l collected 
again as a function of the seniority to explore a possible dependence of x^- We 
can see that, except in the special case of ^4 = 48, which shows a quite large 
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Table 7 

Evolution of (%) 

function of the maximum seniority permitted, s^, in the 
wave functions for all the studied decays. 





^ = 48 


A = 76 


A = 82 


A = 124 


A = 128 


A = 130 


A = 136 





3.2 


1.6 


1.3 


0.5 


0.6 


0.5 


0.4 


4 


9.0 


1.0 


0.8 


0.1 


0.0 


-0.2 


-0.7 


6 


9.4 


0.9 


0.9 


0.1 


0.1 


-0.1 


-0.6 


8 


9.9 


0.6 


0.3 


0.1 


-0.3 


-0.5 


-0.6 


10 




0.6 


0.2 


0.1 


-0.3 


-0.5 




12 




0.4 


0.2 


0.1 


-0.3 






14 




0.4 













ratio, all the other numbers are far from these found in ref. [15], and could be 
considered negligible as in ref. [I6j. This remains true although to a slightly 
lesser extent if we keep only the = 4 figures, which are larger than the 
corresponding complete space values — thus pointing to a similar dependence 
on the seniority truncation to that of the Fermi component — but still very 
minor to be significant in the final result. 

On the other hand, a fully consistent treatment of the SRC's would demand 
regularizing the Oz//3/3 operator using the same prescription than for the bare 
interaction. However, this approach is beyond present ISM or QRPA capabil- 
ities. Hence, general prescriptions, that also come from the regularization of 
bare interactions into the nuclear medium, are used instead. 

The results presented so far have been obtained using the standard Jastrow 
type correlator of eq. ( fTTl) . Other authors [35] have recently argued that this 
correction is somewhat too aggressive, and have proposed another method — 
namely the Unitary Correlation Operator Method (UCOM) [19] — to estimate 
the SRC, which leads to a much smoother correction, of the order of 5% 
compared to the 15-25% of the Jastrow correlator, see Table O 

We have estimated the value of our ISM results taking an UCOM type SRC 
by simulating the correlator as that of the ST = 01 channel [36], common 
throughout the calculation. The correlator of the other important — even — 
channel is very similar to this one, and the difference should not change our 
estimated results. The numbers obtained are listed in Table El where we find 
a 5% reduction of the NME with this UCOM SRC ansatz, in agreement with 
the QRPA calculations. Thus, treating the SRC with a softer prescription of 
this type increases our Jastrow correlated final results by some 20%, leaving 
the full reduction of the NME due to all contributions at around 25-35%. 
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Table 8 

Contributions of the Jastrow and UCOM type SRC to the NME. Jastrow results of 
ref.[l8] were slightly different from present because only one cutoff parameter was 
considered in the FNS terms. 



0^(3/3 Transition 






Jastrow 


AMucom{%) 


Jastrow {%) 






0.92 


0.85 


0.61 


8 


34 






2.96 


2.81 


2.30 


5 


22 




^ 82Kr 


2.79 


2.64 


2.18 


5 


22 


i24Sn- 


. i24Te 


2.77 


2.62 


2.10 


5 


24 


i28Te - 


^ i28xe 


3.05 


2.88 


2.34 


6 


23 


isoTe - 


^ i30xe 


2.81 


2.65 


2.12 


6 


25 


i36Xe - 




2.32 


2.19 


1.76 


6 


24 



In figure [T] the ISM and QRPA results for the NME's are compared within 
this UCOM treatment of the SRC. The same figure but considering Jastrow 
type SRC was shown in ref. [H], but inadvertently the results obtained with 
different values of gA where not rescaled properly. This is corrected in figure 
O where the comparable M'^'^^^'s introduced in eq. ^ are represented. By 
comparing both figures, it is confirmed that there is a common trend; the 
QRPA values are larger than the ISM ones when the high seniority components 
are important in the latter. For both ISM and QRPA the only net effect of 
UCOM is an enhancement of the Jastrow results of about 20%. 

Whether the UCOM or Jastrow method is more appropriate to treat the Oz//?/? 
short range correlations is still an open question. Therefore, taking into ac- 
count the limitations of our method regarding the SRC's, this different results 
obtained by both prescriptions may be considered as an estimation of the 
range of the effect of SRC's —5-20%. 



4 Radial dependence of the NME 

Very recently the radial evolution of the NME has been studied in ref. [20] , 
in order to see for what internucleonic distances r the NME gets the major 
contribution. This is done by representing the operator C(r) defined as: 

oo 

M""^'^ = J C{r)dr (12) 
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Figure 1. The neutrinoless double beta decay M"^''^^'s for ISM and QRPA calcu- 
lations treating the SRC with the UCOM approach. Tu07 QRPA results from ref. 
[20] and Jy07 results from refs. |16jl7j . 



The form of this operator for the transitions studied is shown in Figure [3l 
This result is in full agreement with the QRPA, confirming the findings of ref. 
[20] . This is, beyond r = 3 fm there is no overall contribution to the NME, 
while the maximum value of C{r) occurs around r = 1 fm, which means 
that almost the complete value of the NME comes from the contribution of 
decaying nucleons which are close to each other. This distance corresponds to 
a momentum of g ~ 200 MeV, twice the expected value estimated in Section [21 
Such a small distance is partly due to the cancellation that happens between 
the contribution of decaying pairs coupled to J = and J > 0, as can be seen 
for the dominant GT component of the ®^Se ^^Kr 0i>/3(3 decay in Fig. [4] — for 
all the other transitions the same tendency is reproduced. That the NME's 
radial shape of the QRPA and ISM calculations be identical is quite intriguing 
and perhaps points to a hidden simplicity in their physics. In fact, one can 
argue that QRPA and ISM NME's, radial dependence included, differ only by 
a scaling factor, which can be expressed as the ratio of the average number of 
pairs in both calculations. 
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Figure 2. Same as figlUbut with Jastrow type SRC. Tu07 QRPA results from ref. 
[26] and Jy08 results from ref. |30| . 



5 Decays to Of excited states 



When considering the 0^(3(3 process, the most favourable decay channel is, 
due to phase space considerations, the 0"*" ground state to ground state one, 
which is the case of all transitions so far presented in this paper. However, 
for experimental reasons, it would be very interesting to have a decay to 
an excited state if it is not too much suppressed compared to the decay to 
the , because the background reduction, coming from the simultaneous 
detection of the electrons from the Ou^P decay and the photon(s) from the 
decay of the final nucleus excited state, might make up for the longer lifetime. 
In order to find a candidate for this final excited state, higher spins have to 
be discarded since these decays are disfavored not only by the phase space 
but also by the transition operator, which in order to couple different spin 
states is necessarily of higher order in the current — and hence much smaller 
in magnitude. However, excited Oi" states could have a chance, as they share 
the transition operator with the decays to the ground state. In this new case 
the phase factor would disfavour the decay but, if the NME pushed in the 
opposite direction and increased it enough, then the lifetime of the process 
would not be much larger than the corresponding to the 0^^, . To explore this 
possibility, we have computed, for the first time within the ISM, the NME's 
of the Ou^P decays to excited Oi" states. 
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Figure 3. Evolution of C(r) for the ^^Ca ^ ^^Ti, ^^Se ^^Ki, ^^^Sn ^^^Te, 
^^°Te ^^°Xe and ^^^Xe ^^^Ba transitions. SRC are not included in the calcula- 
tion. 

Table 9 

NME's for the different decays to the different excited 0^ states studied. The ground 
states are denoted by gs. Half-lives are calculated for {niu) = 1 eV. The minor 
differences between these NME's and those of ref. [18] come as in fig. [H 
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The results are gathered in Table [S With the only exception of the relatively 
small increase of our result for A = 48, not really significant because of its 
huge space factor suppression, we see that a common feature of all calculations 
is that the NME's for the decays to O^" states are smaller than the decays to 
the ground state. This has been the case also in previous QRPA calculations. 
However, since these results have been superseeded by more recent ones for 
the initial and final ground states — the case of refs. [2T|22f23p4] . which do 
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Figure 4. Evolution of the GT part of C(r) for the ®^Se ®^Kr transition. The 
contribution of decaying pairs coupled to J = and J > is also shown. SRC are 
not included. 

not consider the HOC relevant terms — or do not treat SRC properly — as 
happens in ref. [25j — , new calculations for the transitions to final 0^ states 
are required to be comparable to our numbers. 

Table M also includes the predicted half-lives for the transitions. We see that 
our results are typically two orders of magnitude longer for the decays to ex- 
cited states. The least disfavored gs — > O^*" transition would be that of ™Ge, 
which is hindered by a factor 2.4 from the NME times 12, the reduction factor 
coming from the phase space. This is, in that case the gs Of transition is 
suppressed by a factor 25 — 30 compared to the gs — > gs, which is probably 
too large to be compensated by the experimental gain via background reduc- 
tion. Nevertheless, it corresponds to experimentalists to evaluate the practical 
interest of the decay to the excited 0"*" in view of such a suppressed rate. 



6 Summary and Conclusions 

In this work we have calculated within the ISM framework the NME for Oz//?/? 
decays of the emitters ^^Ca, ^^Ge, ^^Se, ^^'^Sn, ^^^Te, ^^°Te and ^^exe. Also, we 
have studied the relative importance of different contributions to the NME, 
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namely HOC, FNS and SRC, concluding that all of there is a nice agreement 
between the ISM and different QRPA methods. However, the values of the 
QRPA differ from the ISM ones when the high seniority components are im- 
portant in the latter. We surmise that this difference may be due to the QRPA 
underestimation of the high seniority components of the wave functions in 
these cases. In addition, some discrepancies remain regarding the importance 
of the effect of the Tensor part of the NME, and there is still uncertainties in 
the treatment of the SRC's. 

We have also studied the radial behaviour of the NME's, finding again agree- 
ment between the ISM and the QRPA results. 

Finally, we have calculated the NME's for the decay to excited Of states. 
These transitions are suppressed with respect to the usual ones to the ground 
state to a different extent for each decay, but at least in a factor 25 — 30 for 
the least hindered A — 76. 
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